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Abstract The Lifshitz theory provides a method to calculate the Casimir force be- 
tween two flat plates if the frequency dependent dielectric function of the plates is 
known. In reality any plate is rough and its optical properties are known only to 
some degree. For high precision experiments the plates must be carefully character- 
ized otherwise the experimental result cannot be compared with the theory or with 
other experiments. In this chapter we explain why optical properties of interacting 
materials are important for the Casimir force, how they can be measured, and how 
one can calculate the force using these properties. The surface roughness can be 
characterized, for example, with the atomic force microscope images. We introduce 
the main characteristics of a rough surface that can be extracted from these images, 
and explain how one can use them to calculate the roughness correction to the force. 
At small separations this correction becomes large as our experiments show. Finally 
we discuss the distance upon contact separating two rough surfaces, and explain the 
importance of this parameter for determination of the absolute separation between 
bodies. 
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1 Introduction 

The Casimir force fl\ between two perfectly reflecting metals does not depend on 
the material properties. This is rather rough approximation as the Lifshitz theory 
demonstrates ||2][3]|4| (see also Dzyaloshinskii and Pitaevskii paper in this volume). 
In this theory material dependence of the force enters via the dielectric functions 
of the materials. Because the Casimir-Lifshitz force originates from fluctuations of 
the electromagnetic field, it is related with the absorption in the materials via the 
fluctuation dissipation theorem. The dissipation in the material at a frequency (o is 
proportional to the imaginary part of the dielectric function e{(o) ~ e {(o) + ie (o). 
Thus, to predict the force one has to know the dielectric properties of the materials. 

In most of the experiments where the Casimir force was measured (see reviews 
||5] |6l and Lamoreaux paper in this volume) the bodies were covered with conduct- 
ing films but the optical properties of these films have never been measured. It is 
commonly accepted that these properties can be taken from tabulated data in hand- 
books QUI. Moreover, for conducting materials one has to know also the Drude 
parameters, which are necessary to extrapolate the data to low frequencies ||9l. This 
might still be a possible way to estimate the force, but it is unacceptable for calcu- 
lations with controlled precision. The reason is very simple ifTOl [TTl [T2HT31 : optical 
properties of deposited films depend on the method of preparation, and can differ 
substantially from sample to sample. 

Analysis of existing optical data for Au lfT4l revealed appreciable variation of the 
force in dependence on the optical data used for calculations. Additionally, we mea- 
sured our gold films using ellipsomety in a wide range of wavelengths 0.14 — 33/im 
ifTSll . and found significant variation of optical properties for samples prepared at 
different conditions. Considerable dependence of the force on the precise dielectric 
functions of the involved materials was also stressed for the system solid-liquid- 
soUd JT6I . 

The Lifshitz formula can be applied to two parallel plates separated by a gap d. 
In reality each plate is rough and the formula cannot be applied directly. When the 
separation of the plates is much larger than their root-mean-square (rms) roughness 
w one can calculate correction to the force due to roughness using the perturba- 
tion theory. But even in this case the problem is far from trivial ifTTl [Tsl [T9l . The 
roughness correction can be easily calculated only if one can apply Proximity Force 
Approximation (PFA) 1201 . Application of this approximation to the surface pro- 
file is justified when this profile changes slowly in comparison with the distance 
between bodies. Typical lateral size of a rough body is given by the correlation 
length ^ . Then the condition of applicability of PFA is ^ ^ d. This is very restric- 
tive condition since, for example, for thermally evaporated metallic films the typical 
correlation length is ^ 50 nm. 

The roughness of the interacting bodies restricts the minimal separation do be- 
tween the bodies. This distance (distance upon contact) has a special significance 
for adhesion, which under dry conditions is mainly due to Casimir/van der Waals 
forces across an extensive noncontact area 1211 . It is important for micro (nano) 
electro mechanical systems (MEMS) because stiction due to adhesion is the ma- 



Characterization of optical properties 



3 



jor failure mode in MEMS Il22l . Furthermore, the distance upon contact plays an 
important role in contact mechanics, is very significant for heat transfer, contact 
resistivity, lubrication, and sealing. 

Naively one could estimate this distance as the sum of the rms roughnesses of 
body 1 and body 2, do ^ wi + W2 ll23l . however, the actual minimal separation is 
considerably larger. This is because do is determined by the highest asperities rather 
then those with the rms height. An empirical rule 1241 for gold films gives do ~ 
3.7 X (wi + W2) for the contact area of 1 jim^. The actual value of do is a function 
of the size of the contact area L. This is because for larger area the probability to 
find a very high peak on the surface is larger 

Scale dependence (dependence on the size L) is also important for the Casimir 
force in the noncontact regime. In this case there is an uncertainty in the separation 
5d{L), which depends on the scale L. The reason for this uncertainty is the local 
variation of the zero levels, which define the mathematical (average) surfaces of 
the bodies. This uncertainty depends on the roughness of interacting bodies and 
disappears in the limit L — > 00. 

In this paper we explain how one can collect the information about optical prop- 
erties of the materials, which is necessary for the calculation of the Casimir-Lifshitz 
force. It is also discussed how the optical spectra of different materials manifest 
themselves in the force. We introduce the main characteristics of rough surfaces and 
discuss how they are related with the calculation of the roughness correction to the 
force. Scale dependence of the distance upon contact is discussed, and we explain 
significance of this dependence for the precise measurements of the force. 



2 Optical properties of materials and the Casimir force 

Most Casimir force measurements were performed between metals Il5ll6l l25ll26ll27l 
either e-beam evaporated or plasma sputtered on substrates. For such metallic films 
the grains are rather small in the order of tens of nanometers, and the amount of 
defects and voids is large ll28l . The force measured between silicon single crystal 
and gold coated sphere ||29l simplify situation only partly: the optical properties of 
the Si-crystal are well defined but properties of Au coating are not known well. 

A detailed literature survey performed by Pirozhenko et al. 11141 revealed signif- 
icant scatter in the dielectric data of gold films collected by different groups. The 
measurement errors were not large and could not explain the data scattering. It was 
concluded that scattering of the data for gold films could lead to uncertainty in the 
calculated force up to 8% at separations around 100 nm. Most of the optical data for 
metals do not extend beyond the wavelength of 14 /im in the infrared range ll30l[3TI . 
Thus, it would be important to explore more the infrared regime and compare mod- 
ern measured optical properties of samples used in force measurements with the old 
data. Moreover, mid and far infrared data are very important for the force prediction 
(see Sec. l2.2.2] i. This was accomplished by Svetovoy et al. [151 where ellipsometry 
from the far infrared (IR) to near ultraviolet (UV) was used over the wavelength 
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range 140 nm — 33 /im to obtain the frequency dependent dielectric functions for 
gold films prepared in different conditions. Analysis of different literature sources 
where the dielectric functions of a number of dielectrics such as silica and some 
liquids was performed by van Zwol et al. lITSl . Situations where the data scattering 
can change even the qualitative behavior of the force (from attractive to repulsive) 
were indicated. 



2.1 Dielectric function in the Casimir force 



In this section we discuss how the dielectric functions of materials enter the Lifshitz 
theory and how these functions can be found experimentally. 



2.1,1 The force 

Let us start the discussion from the Lifshitz formula Q between two parallel plates 
separated by a gap d. It has the following form 

. , kT ^ , 7 „ rJ'ri'e-^^co^/ 

where "prime" at the sign of sum means that the n = term must be taken with the 
weight 1/2, the wave vector in the gap is K = (q, Kq) with the z-component fCo de- 
fined below. Index "0" is related with the gap. Here rj'j are the reflection coefficients 
of the inner surfaces of the plates (index 1 or 2) for two different polarizations: v = s 
or transverse electric (TE) polarization, and v = p ot transverse magnetic (TM) po- 
larization. Specific of the Lifshitz formula in the form ^ is that it is defined for a 
discrete set of imaginary frequencies called the Matsubara frequencies 

iKkT 

(On^iCj^i— — «, « = 0,1,2,..., (2) 
« 

where T is the temperature of the system and n is the summation index. 

In practice the interacting bodies are some substrates covered with one or a few 
layers of working materials. If the top layer can be considered as a bulk layer then 
the reflection coefficients for body / are given by simple Fresnel formulas 1321 : 

^ ^ Kp-Kj ^ eijiQKQ-eQjiQKi 
' Kb + K-,' ' eiiiQKQ + eoiiQKi' 

where 



Ko = \jeo{iQ^+q^, Ki = Jei{iQ^+q^. (4) 
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For multilayered bodies these formulas can be easily generalized (in relation with 
the dispersive forces see Ref. |f33]| ). Only the reflection coefficients depend on the 
dielectric functions of the plate materials; the function eo('C) of the gap material 
enters additionally in Kq. 

At small separations the thermal dependence of the force is very weak and in 
many cases can be neglected. Because important imaginary frequencies are around 
the so called characteristic frequency = c/2d, then the relative thermal correction 
can be estimated as kT /h^c- For room temperature T ~ 300° K and separations 
smaller than 100 nm the correction will be smaller than 3%. If one can neglect this 
correction then in the Lifshitz formula ^ can be considered as a continuous variable 
and the sum in ^ is changed by the integral according to the rule: 



The material function e(/Q (we suppress the indexes for a while) cannot be 
measured directly but can be expressed via the observable function £ (o) with the 
help of the Kramers-Kronig (KK) relation ll32l : 



The knowledge of e(/Q is of critical importance for the force calculations. Equation 
^ demonstrates the main practical problem. To find the function e(/Q for ^ ~ 
in general one has to know the physical function £ (co) in a wide range of real 
frequencies, which is not necessary coincides with co ^ ^c- This property of the 
Casimir force was stressed in Ref. iT2il and then was demonstrated experimentally 
Il34ll35ll36l . It win be discussed below on specific examples. 

2.1.2 The optical data 

The optical properties of materials are described by two measurable quantities: the 
index of refraction n{X) and the extinction coefficient k{X), which both depend on 
the wavelength of the electromagnetic field X . Combined together they define the 
complex index of refraction n{X) = n{X) + ik{X). The real part defines the phase 
velocity in a medium v ~c/n where c is the speed of light. The imaginary part tells 
us how much light is adsorbed when it travels through the medium. The dielectric 
response of a material for the UV (/io) > 5 eV), IR (0.01 — 1 eV) and MicroWave 
(MW) or TeraHertz range (10^^ — 10^^ eV), is related to electronic polarization 
resonances, atomic polarization resonances (in case of metals this is a gas of quasi 
free electrons), and dipole relaxation, respectively. 

The complex dielectric function e(fi)) = £■'{(£>) + i£"{(£>) and the complex index 
of refraction are related as e(fi)) = «'(fi)), which is equivalent to the following equa- 




(5) 




(6) 
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tions: 

e'=n^-k^, e" = 2nk. (7) 

In many cases only the absorbance is measured for a given material. In this case the 
refraction index can be found from the KK relation at real frequencies: 

e{(o)^\ + -P ( dx^^^, (8) 



where P means the principal part of the integral. 

Kramers-Kronig relations originating from causality have a very general char- 
acter. They are useful in dealing with experimental data, but one should be careful 
since in most cases dielectric data are available over limited frequency intervals. As 
a result specific assumptions must be made about the form of the dielectric data out- 
side of measurement intervals, or the data should be combined with other (tabulated) 
experimental data before performing the KK integrals. 

A powerful method to collect optical data simultaneously for both e and e is el- 
lipsomery. This is a non destructive technique where one measures an intensity ratio 
between incoming and reflected light and the change of the polarization state. Ellip- 
sometry is less affected by intensity instabilities of the light source or atmospheric 
absorption. Because the ratio is measured no reference measurement is necessary. 
Another advantage is that both real and imaginary parts of the dielectric function 
can be extracted without the necessity to perform a Kramers-Kronig analysis. The 
ellipsometry measures two parameters f and /i, which can be related to the ratio of 
complex Fresnel reflection coefficients for p- and s-polarized light 13711381 

p = — =\.imWe'^, (9) 

where r^'* are the reflection coefficients of the investigated surface, and the angles 
f and A are the raw data collected in a measurement as functions of the wave- 
length X. When the films are completely opaque (bulk material), then the reflection 
coefficients are related with the dielectric function as follows 



(e) cos & - \ (e) - sin cos \ (e) — sin 

r„ = \ = / (10) 

(e) cos + W (e) - sin cos + W (e) — sin i> 

where t? is the angle of incidence and (e) = (e(A)) is the "pseudo" dielectric func- 
tion of the films. The term "pseudo" is used here since the films may not be com- 
pletely isotropic or uniform; they are rough, and may contain absorbed layers of 
different origin because they have been exposed to air. If this is the case then the di- 
electric function extracted from the raw data will be influenced by these factors. The 
dielectric function is connected with the ellipsometric parameter p for an isotropic 
and uniform solid as 
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( 



1-p 



) 



i 



e = sin^ 1} 1 + tan^ 



1+P 



(11) 



As it was stated before the spectral range of our measured data is from 137 nm 
to 33 jUm. Even longer wavelengths have to be explored to predict the force be- 
tween metals without using the extrapolation. Ellipsometry in the terahertz range 
0.1—8 THz (wavelengths 38 — 3000 /im) is difficult due to lack of intense sources 
in that range, and these systems are still in development 1391 . Typically synchrotron 
radiation is used as a source deeming these measurements very expensive. Nonethe- 
less for gold films it would be extremely interesting to have dielectric data in this 
regime. 

Dielectric data obtained by ellipsometry or absorption measurements f40ll in the 
far UV regime are also rare. The most obvious reason for this is that these mea- 
surements are expensive because high energy photons must be produced, again at 
synchrotrons BTI . Furthermore, ellipsometry in this range becomes complicated as 
polarizing materials become non transparent. For this range a few ellipsometry se- 
tups exist around the world covering the range 5 — 90 eV (wavelengths 12 — 200 nm) 
ll4ll . The vacuum UV (VUV) and extreme UV (XUV) parts may not be very im- 
portant for metals but for low permittivity dielectrics such as all liquids, and, for 
example, silica or teflon, there is a major absorption band in the range 5 — 100 eV 
(see Fig.|6l). It is precisely this band that dominates in the calculations of the Casimir 
force for these materials. It is very unfortunate that precisely for this band dielec- 
tric data are lacking for most substances except for a few well know cases as, for 
example, water. 



In this section we discuss optical characterization of our gold films prepared in dif- 
ferent conditions using ellipsometers. Then we discuss the dielectric function at 
imaginary frequencies for Au films and for metals in general stressing the impor- 
tance of very low real frequencies for precise evaluation of £{iQ at i^^ ^ 1 eV 
(separations around 100 nm). Finally we describe variation in the Casimir force if 
different samples would be used for the force measurements. 

2.2.1 £ (to) for Au films 

Let us have now a closer look at the dielectric functions of our gold films ifTSll used 
for the force measurements in Refs. Il42ll24l . For optical characterization we have 
prepared five films by electron beam evaporation. Three of these films of different 
thicknesses 100, 200 and 400 nm were prepared within the same evaporation system 
on Si with 10 nm titanium sublayer and were not annealed. Different evaporation 
system was used to prepare two other films. These films were 1 20 nm thick. One film 
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Fig. 1 Flattened roughness scans (up to 4000 x 4000 pixels) of gold surfaces where the highlighted 
areas are magnified. The color scale bars can be applied only to the large images, (a) 100 nm Au 
on Si, (b) Au coated polysterene sphere (first plasma sputtered then 100 nm Au evaporated, (c) 
1600 nm Au on Si, (d) very high quality 120 nm Au on mica, annealed for a few hours and slowly 
cooled down. 

was deposited on mica and was extensively annealed. The other one was deposited 
on Si with chromium sublayer and was not annealed. 

The AFM scans of the 100 nm film and the annealed 120 nm film on mica are 
shown in Fig. [T] In the same figure are shown also the gold covered sphere and 
1600 nm film, which where not used in optical characterization. One can see that 
the annealed sample is atomically smooth over various length scales with atomic 
steps and terraces visible. Nevertheless, the local trenches of 5 nm deep are still 
present. 

Optical characterization of the films was performed by J. A. Woollam Co., Inc. 
ll43l . Vacuum ultraviolet variable angle spectroscopic ellipsometer (VASE) was used 
in the spectral range 137 — 1698 nm. In the spectral range 1 .9 — 32.8 jJ-m the infrared 
variable angle spectroscopic ellipsometer (IR-VASE) was used for two incidence 
angle of 65° and 75°. The roughness and possible absorbed layer on the film sur- 
face can have some significance in the visible and ultraviolet ranges but not in the 
infrared, where the absorption on free electrons of metals is very large. Moreover, 
the effect of roughness is expected to be small since for all films the rms roughness 
is much smaller than the smallest wavelength 137 nm. Because the infrared domain 
is the most important for the Casimir force between metals, we will consider (e(A)) 
extracted from the raw data as a good approximation for the dielectric function of a 
given gold film. 
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Fig. 2 (a) Measured e" as a function of frequency ft), (b) The same for | e' | . For clearness the results 
are presented only for three samples 2, 3, and 5. 

Figure|2a) shows the experimental results for e" (o) for three of five investigated 
samples. Around the interband transition (minimum of the the curves) the smallest 
absorption is observed for the sample 5 (annealed on mica) indicating the smallest 
number of defects in this sample 1281 . On the contrary, this sample shows the largest 
I e (o) I in the infrared as one can see in Fig. [2}5. An important conclusion that can be 
drawn from our measurements is the sample dependence of the dielectric function. 
The sample dependence can be partly attributed to different volumes of voids in 
films as was proposed by Aspnes et al. Il28l . The values of e and their dispersion for 
different samples are in good correspondence with old measurements ll30l[3TI . The 
log-log scale is not very convenient for having an impression of this dependence. 
We present in Table I the values of e for all five samples at chosen wavelengths 
A = 1,5, 10 jUm. One can see that the real part of £ varies very significantly from 
sample to sample. 



Sample 


X = 1 /im 


A = 5 tim 


A = 10 ^m 


l,400nm/Si 


-29.7 + /2.1 


-805.9 + ;185.4 


-2605.1 +(T096.3 


2, 200 nm/Si 


-31.9 + ;2.3 


-855.9 + ;195.8 


-2778.6 + (T212.0 


3, 100 nm/Si 


-39.1 + ;2.9 


-1025.2 + (264.8 


-3349.0 + (T574.8 


4, 120 nm/Si 


-43.8 + ;2.6 


-1166.9 + (213.9 


-3957.2 + (1500.1 


5, 120 nm/mica 


-40.7 + /1. 7 


-1 120.2 + (T78.1 


-4085.4 + (1440.3 



Table 1 Dielectric function for different samples at fixed wavelengths A = 1,5, 10^ m. 



One could object that the real part of £ does not play role for £(/C) and only vari- 
ation of e' (ft)) from sample to sample is important. However, both £ ' (co) and £ (o) 
are important for precise determination of the Drude parameters. Let us discuss now 
how one can extract these parameters from the data. 

All metals have finite conductivity. It means that at low frequencies (O ^ Q the 
dielectric function behaves as £{(£)) -^Ako/o), where a is the material conductivity. 
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It has to be stressed that this behavior is a direct consequence of the Ohm's law and, 
therefore, it has a fundamental character Because the dielectric function has a pole 
at 0) — > 0, the low frequencies will give a considerable contribution to e{iQ even if 
^ is high (for example, in visible part of the spectrum) as one can see from Eq.®. 

Usually it is assumed that at low frequencies the dielectric functions of good 
metals follow the Drude model: 

e(fi)) = l , ^ . (12) 

where cOp is the plasma frequency and 7 is the relaxation frequency of a given metal. 
When (B — !> we reproduce the 1/(0 behavior with the conductivity a = (Op/A-ny- 
For good metals such as Au, Ag, Cu, Al typical values of the Drude parameters are 
(Op ~ 10'^ rad/s and 7- lO'"* rad/s. 

Separating real and imaginary parts in Eq. (fT2b one finds for e and e 

/ (ol „ (oly 

e(a)) = l- £ (0))= . / . (13) 

These equations can be applied below the interband transition (O < 2.45 eV {X > 
0.5 jim) |44], but this transition is not sharp and one has to do analysis at lower 
frequencies. Practically Eq. ( fT3] l can be applied at wavelengths X > 2 /im that co- 
incides with the range of the infrared ellipsometer The simplest way to find the 
Drude parameters is to fit the experimental data for e and e with both equations 
( fT3] ). Alternatively to find the Drude parameters one can use the functions n{(o} and 
k{co) and their Drude behavior, which follows from the relation e(a)) =n^(a)). This 
approach uses actually the same data but weight noise differently. 

Completely different but more complicated approach is based on the KK rela- 
tions (O (see Refs. |[T4| and ffTSl for details). In this case one uses measured e (©) 
extrapolated to low frequencies according to the Drude model and extrapolated to 
high frequencies as A/(o^, where A is a constant. In this way we will get e (o) 
at all frequencies. Using then Eq. dH) we can predict e (to). Comparing the predic- 
tion with the measured function we can determine the Drude parameters. Similar 
procedure can be done for n{(o) and k{(o). 





1, 400 nra/Si 


2, 200 nm/Si 


3, 100 nm/Si 


4, 120 nm/Si 


5, 120 nm/mica 


7 (meV) 


40.5 ±2.1 


49.5 ±4.4 


49.0 ±2.1 


35.7±5.1 


37.1 ± 1.9 


m„ (eV) 


6.82 ±0.08 


6.83 ±0.15 


7.84 ±0.07 


8.00±0.16 


8.38±0.08 


£, (nm) 


22 


26 


32 


70 


200 


w (nm) 


4.7 


2.6 


1.5 


1.5 


0.8 



Table 2 Drude parameters 7, cOp and roughness parameters, the correlation length ^ and rms 
roughness w, for all five measured samples. 
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Fig. 3 The infrared data as functions of the wavelength A for e' and e " (solid lines) and the best 
Drade fits (dashed lines) for two gold films. Panel (a) shows the data for annealed sample 5 and 
panel (b) shows the same for unannealed sample 3. 



All methods for determination of the Drude parameters give reasonably close 
values of both parameters. We cannot give preference to any specific method. In- 
stead, we average the values of the parameters determined by different methods, 
and define the rms error of this averaging as uncertainty in the parameter value. The 
averaged parameters and rms errors are given in the Table |2] We included in this 
table also the correlation lengths £, and the rms roughness w for the sample rough- 
ness profiles. One can see the COp and ^ correlate with each otherQ. This correlation 
has sense because ^ describes the average size of the crystallites in the film; the 
larger the crystallites the smaller number of the defects has the film and, therefore, 
the larger value of the plasma frequency is realized. 

Quality of the Drude fit one can see in Fig.[3]for samples 3 and 5. The fit is prac- 
tically perfect for high quality sample 5, but there are some deviations for sample 3 
at short wavelengths especially visible for e . More detailed analysis ifTSl revealed 
presence of a broad absorption peak of unknown nature around A = 10 /zm. The 
magnitude of this absorption is the largest for poor quality samples 1 and 2. 

2.2.2 e{iQ for metals 

Any method of optical characterization has a minimal accessible frequency cOait 
(cut-off frequency). At co > (Ocut one can measure the dielectric function but at lower 
frequencies (o < (Ocut ont has to make an assumption on the behavior of £(o), i. e. 
extrapolate to low frequencies. In KK relation (|6]l one can separate two intervals: 
(0 < (Ocut, where e (co) has to be extrapolated and (O > (Ocut, where e (o)) is mea- 
sured. Then we can present e(/Q as 



' This congelation was not noted in 1151 and stressed here for the first time. 
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Fig. 4 (a) Relative contribution e„„[iQ/£{il^) to the dielectric function of gold at imaginary fre- 
quencies originating from the extrapolated region CO < COcut (see explanations in the text), (b) The 
dielectric functions at imaginary frequencies for samples 2,3, and 5; the thick curve marked as 
corresponds to an "ideal sample" with the plasma frequency of single crystal. 



2 y (oe {a) 2 / coe o) 

EcutilQ^- / d(0 ^ ' eexper{lQ = - / d(0 ' ' . (14) 
7Z J CO^ + Q^ K J CO^ + C 

(0„„ 

Of course, at very high frequencies we also do not know e" (o) but, for metals high 
frequencies are not very important. For this reason we include the high frequency 
contribution to SexperiiQ- 

We can estimate now the contribution of EcutiiQ to £('C)- For that we assume 
the Drude behavior at o < (Oa,t with the parameters cOp ~ 9.0 eV and 7 = 35 meV 
At higher frequencies co > (Ocm we take the data from the handbook Q, for 
which the cut-off frequency is fi)t,„ = 0.125 eV. These extrapolation and data were 
used for interpretation most of the experiments, where the Casimir force was mea- 
sured. In Fig. |4{a) the solid curve is the ratio EcutiiQ/siiQ calculated with these 
data. One can see that at ^ = 1 eV (c/ ^ 100 nm) the contribution to £{iQ from 
the frequency range (O < (Oait is 75%. It means, for example, that if we change the 
Drude parameters, three fourths of e{iC) will be sensitive to this change and only 
one forth will be defined by the measured optical data. Therefore, the extrapolation 
procedure becomes very important for reliable prediction of £{iQ- 

The Drude parameters can vary from sample to sample due to different density 
of defects. The plasma frequency is related with the density of quasifree electrons 

as (Op ~ 4-nNe^/m*, where for gold m* ss me is the effective mass of electron. 
The value (Op = 9.0 eV is the maximal value of this parameter, which corresponds 
to in a single crystal Au. In this way (Op was estimated in Ref. All deposited 
films have smaller values of cOp as one can see from Tab. |2] because the density of 
the films is smaller than that for the single crystal material. Precise values of the 
Drude parameters are extremely important for evaluation of £(;Q and finally for 
calculation of the force. 

The dielectric function e(/Q becomes less dependent on the Drude parameters 
if the cut-off frequency is smaller. For example, our optical data 115 1 were collected 
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up to minimal frequency cOcut = 38 meV that is about four times smaller than in 
the handbook data. The dashed curve in Fig.|4ja) shows the ratio Scut/s calculated 
for our sample 3 with the Drude parameters cOp = 7.84 eV and 7 = 49 meV. Now 
the relative contribution of EcutiiQ at ^ = 1 eV is 37%. It is much smaller than for 
handbook data, but still dependence on the precise Drude parameters is important. 
Let us imaging now that we have been able to measure the dielectric response of 
the material for the same sample 3 from iflSl to frequencies as low as 1 THz. In 
this case the cut-off frequency is = 4 meV and the relative contribution of the 
extrapolated region £cut/£ is shown in Fig.Ua) by the dash-dotted line. Now this 
contribution is only 5% at 1^ = 1 eV. 

The dielectric functions £i{iQ, where / = 1,2, ..,5 is the number of the sample, 
were calculated using the Drude parameters from Tab. |2] As a reference curve we 
use eo('C)' which was evaluated with the parameters cOp ~ 9.0 eV and 7= 35 meV 
in the range o < 0.125 eV and at higher frequencies the handbook data Q were 
used. The results are shown in Fig. |5b). As was expected the maximal dielectric 
function is eo('C)' which corresponds in the Drude range to a perfect single crystal. 
Even for the best sample 5 (annealed film on mica) the dielectric function is 15% 
smaller than eo('C) at 1^ = 1 eV. For samples 1 and 2 the deviations are as large as 
40%. 



2.2.3 The force between Au films 

It is convenient to calculate not the force itself but so called the reduction factor 77, 
which is defined as the ratio of the force to the Casimir force between ideal metals: 

r^id) = -^, F^id) = ~^. (15) 
' F':{dy ^ ' 240c/4 

We calculate the force between similar materials at T = using the substitute (|5]i in 
Eq. ([1]). For convenience of the numerical procedure one can make an appropriate 
change of variables so that the reduction factor can be presented in the form 

where the reflection coefficients as functions of x and y are defined as 



1 -s _ £{i^cxy)-s 



with 

c 



= ^l+x^£iiCexy)-l], Cc = ^- (18) 
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Fig. 5 (a) Reduction factor T] as a function of tiie separation d for samples 1, 3, and 5. The thick line 
shows the reference result calculated with eo(;C)- (b) Relative deviations of the reduction factors 
for different samples from the reference curve T]()(rf), which were evaluated using the handbook 
optical data Q and the Drude parameters (Op = 9 eV, tt)j = 35 meV. 



The integral ( fT6] l was calculated numerically with different dielectric functions 
e,(/Q. The results are presented in Fig.|5la) for samples 1, 3, and 5. The reference 
curve (thick line) calculated with eo('C) is also shown for comparison. It repre- 
sents the reduction factor, which is typically used in the precise calculations of the 
Casimir force between gold surfaces. One can see that there is significant difference 
between this reference curve and those that correspond to actual gold films. To see 
the magnitude of the deviations from the reference curve, we plot in Fig. |5jb) the 
ratio (t7o ^ J];) /tjo as a function of distance d for all five samples. 

At small distances the deviations are more sensitive to the value of cOp. At large 
distances the sample dependence becomes weaker and more sensitive to the value 
of (Or- For samples 1 and 2, which correspond to the 400 nm and 200 nm films 
deposited on Si substrates, the deviations are especially large. They are 12-14% at 
J < 100 nm and stay considerable even for the distances as large as 1 fim. Samples 3, 
4, and 5 have smaller deviations from the reference case but even for these samples 
the deviations are as large as 5-7%. 



2.3 Low dielectric materials 



The Lifshitz theory predicts ||3] that the dispersive force can be changed from at- 
tractive to repulsive by choosing the interacting materials immersed in a liquid. Re- 
cently this prediction was confirmed experimentally 1451 (see Capasso paper in this 
volume). Repulsive forces arise when the dielectric functions at imaginary frequen- 
cies in the liquid gap, eo{iQ, is in between the functions of the interacting bodies 1 
and 2: Ei {iQ > eo('C) > £2('C)- One can expect significant dependence on precise 
dielectric functions nearby the transition from attractive to repulsive force. This sit- 
uation is exactly the case for the system silica-liquid-gold. In this section we present 
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Fig. 6 (a) Dielectric data of the materials obtained from references in text, (b) Dielectric functions 
at imaginary frequencies. The solid and dashed lines for silica and gold are two different sets of 
optical data. For water the solid line is from the data in Ref. 1491 . and the dashed line is an 1 1 -order 
oscillator model 1501 . which has been fitted to a different set of optical data. 



calculations for multiple liquids with various degrees of knowledge of the dielectric 
functions. 

Liquids do not have grains or defects, but the density of a liquid is a function of 
temperature 1461 . and as a result the number of absorbers in the liquid varies with 
temperature. Furthermore, liquids can contain impurities like salt ions which can 
change the dielectric function (see discussion in Ref. |47l). Although for metals the 
dielectric function is very large in the IR regime, for liquids and glasses it is not 
the case. Consequently for low dielectric materials the UV and VUV dielectric data 
have a strongest effect on the forces. 

For gold, silica, and water the dielectric functions are well known and measured 
by various groups. Let us consider first the interaction in the system gold-water- 
silica. We will use two sets of data for gold from the previous subsection (sample 1 
and the "ideal sample"). Also two sets of data will be used for silica as obtained by 
different groups BSl . Finally, for water we will use the data of Segelstein compiled 
from different sources 1491 . and an 11 -order oscillator model |f50l that has been fit 
to different sets of data IST] |52l . All the dielectric data are collected in Fig. |6la). 
The corresponding functions at imaginary frequencies are shown in Fig. |6jb). One 
can see considerable difference between solid and dashed curves corresponding to 
different sets of the data. 

It has to be noted that e {(o) for water and silica are very close in a wide range 
of frequencies 5 ■ 10^^ < O < 5 • 10^ eV. As the result at imaginary frequencies 
8si02 ('O ™d £//2o('C) differ on 30% or less in the range 10^^ < C < 10^ eV, which 
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Fig. 7 (a) Variation of the force relative to Fo{d) in the system gold-water-silica. Circles, squares, 
triangles and stars mark the curves which were calculated using different sets of the dielectric data, 
(b) Variation of the force for silica-water-silica (dashed line) with different sets of data for SiOa 
and for gold-water-gold (solid line) using different data for Au. 



is comparable with the magnitude of variation of e(/'Q due to data scattering. This 
similarity in the dielectric functions results in a strong dependence of the Casimir 
force in the system gold-water-silica on the used optical data. It is illustrated in 
Fig. |7ta) where the relative change of the force is shown. The scatter of the force 
riches of the level of 60% for separation d < 500 nm. Even more clear the effect can 
be seen in Fig. |7jb). The solid curve shows variation of the force in Au-water-Au 
system when different optical data for Au are used. In this case the relative change 
of the force is not very large. However, for the system silica-water-silica the use of 
different optical data for silica influence the force very significantly (dashed curve). 

We have to conclude that comparison of force measurements with prediction of 
the Lifshitz theory becomes reliable when the dielectric properties of the specific 
samples used in force measurement are measured over a wide range of frequencies. 

In most of the papers where liquid gap between bodies is studied the dielectric 
function of the liquid is approximated using the oscillator models Il53ll54l . For il- 
lustration purposes we mention that alcohols (and other liquid substances) can be 
described, for example, by a three oscillator model for the dielectric function e{iQ 
at imaginary frequencies 1531 

eiiQ = 1 + -^^^-^ + ^ + ^ (19) 

1 + C/«MW 1 + (C/«/k)2 l + (C/£0t/v)2 

Here no is the refractive index in the visible range, Eq is the static dielectric constant, 
and EjR is the dielectric constant where MW relaxation ends and IR begins. The 
parameters (Omw, (Oir, and (Ouv are the characteristic frequencies of MW, IR, and 
UV absorption, respectively. It has to be stressed that oscillator models should be 
used with caution, because some of them are of poor quality ||T6l . 

For ethanol rather detailed information on the dielectric function exists, but even 
in this case variation in dielectric data was found llT6l l40l . An interesting fact for 
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Fig. 8 Dielectric data at real (a) and at imaginary frequencies (b) for methanol, ethanol, propanol, 
and butanol, as described in the text. For comparison the dielectric data for silica and gold samples 
are also shown in (b). In (c) we show the forces for the gold-alcohol-silica, for two different sets 
of measured dielectric data for silica. In the case of butanol the force is attractive for one set, and 
repulsive for another one. 



higher alcohols is that the absorption in the UV range increases when increasing the 
aUcane chain. In Fig. [8] we show the dielectric data for the first four alcohols. 

The VUV data were taken from ||55l . These measurements were done in the gas 
phase, but they can be converted to the liquid case by considering the number of 
absorbers in gas and liquid. The near UV data was taken from 1561 . For the XUV 
we have only data for ethanol and propanol ISTll . For methanol and butanol we used 
cubic extrapolation, £ ^ 1 / (O^, which is in very good agreement with the cases of 
ethanol and propanol. In the near IR (NIR) to visible (VIS) ranges the extinction 
coefficient k of ethanol, and other alcohols, is very low and can be taken to be zero, 
k^Q, which is qualitatively consistent with the fact that all alcohols are transparent 
in the visible range. The IR data can be found in Ref. 1581 . The far IR (FIR) data are 
known only for methanol 1591 . For the other alcohols we take the similar functional 
behavior as for methanol but with different parameters and extrapolate the data to 
far IR in this way. 

If one has to estimate the dielectric functions for some alcohols, first of all one 
has to have measured data in the range of major IR peaks and even more impor- 
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tantly the UV absorption must be carefully measured. Thus the dielectric functions 
at imaginary frequencies should be reasonably accurate to within the scatter of the 
data as it was found for ethanol 1401 . 

With the optical data for alcohols we calculated the forces in the system gold- 
alcohol-silica. The forces are attractive and become weaker for methanol, ethanol 
and propanol. For butanol they are extremely weak, but still either repulsive or at- 
tractive. Caution is required in the analysis of optical properties in liquids since 
in general the KK consistency has to be applied properly in order to correct for 
variation of the dielectric properties observed in between different measuring se- 
tups. Effectively the force for gold-butanol-silica is screened to within the scatter 
of the forces related to sample dependence of the optical properties of silica. Mea- 
surements between gold and glass with simple alcohols were performed, but experi- 
mental uncertainty, and double layer forces prevented the measurement of this effect 



3 Influence of surface roughness on the Casimir-Lifshitz force 

The Lifshitz formula ([T]l does not take into account inevitable roughness of the in- 
teracting bodies. When rms roughness of the bodies is much smaller than the sep- 
aration, then the roughness influence on the force can be calculated using the per- 
turbation theory iflTl ITSl IT9l . However, when the separation becomes comparable 
with the roughness the perturbation theory cannot be applied. It was demonstrated 
experimentally |[24l that in this regime the force deviates significantly from any the- 
oretical prediction. The problem of short separations between rough bodies is one 
of the unresolved problems. In this section we give introduction into interaction of 
two rough plates and a sphere and a plate. 



3.1 Main characteristics of a rough surface 

Suppose there is a rough plate which surface profile can be described by the func- 
tion li{x,y), where x and y are the in-plane coordinates. An approximation for this 
function provides, for example, an AFM scan of the surface. It gives the height lijj 
at the pixel position xi = A ■ i and yj — A ■ j, where i,j ~ l,2,...,N and A is the 
pixel size related with the scan size as L = A - N. We can define the mean plane of 
the rough plate as the averaged value of the function h{x,y): Ji = A^' / dxdyh{x,y), 
where A is the area of the plate. This definition assumes that the plate is infinite. In 
reality we have to deal with a scan of finite size, for which the mean plane is at 



(601. 




(20) 
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The difference h — /!„,,, although small, is not zero and is a random function of the 
scan position on the plate. This difference becomes larger the smaller the scan size 
is. Keeping in mind this point, which can be important in some situations (see be- 
low), we can consider ( l20b as an approximate definition of the mean plane position. 

An important characteristic of the a rough surface is the rms roughness w, which 
is given as 

'^ = ];^'L i^'^i^i^yj) - • (21) 

It has the meaning of the surface width. More detailed information on the rough 
surface can be extracted from the height-difference correlation function defined for 
the infinite surface as 

8{R) = jj dxdy [/i(r + R)-/i(r)]2, (22) 

where r = [x^y) and R = r' — r. 

A wide variety of surfaces, as for example, deposited thin films far from equi- 
librium, exhibit the so called self-affine roughness which is characterized besides 
the rms roughness amplitude w by the lateral correlation length ^ (indicating the 
average lateral feature size), and the roughness exponent < // < 1 lIM] |62l |63l . 
Small values of // ~ corresponds to jagged surfaces, while large values of H ^ \ 
to a smooth hill vaUey morphology. For the self-affine rough surfaces g{R) scales as 

«(«)-{£";">>^: 

The parameters w, ^ and H can be determined from the measured height-difference 
correlation function g{R). This function can be extracted approximately from the 
AFM scans of the surface. 

To find roughness correction to the force one has to know (see below) the spectral 
density (j{k) of the height-height correlation function C{R). The latter is related with 
g{R) as g{R) ~ 2w^ — C{R). An analytic form of the spectral density for a self-affine 
surface is given by 1641 



^W = — T7^' C= . ... (24) 



Here C is a normalization constant ll64ll63]| and kc = 2n/Lc is the cutoff wavenum- 
ber 



3.2 Roughness correction 



While the separation between two surfaces is large in comparison with the rms 
roughness, d ':$> w, one can use the perturbation approach to calculate the rough- 
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ness correction to the Casimir force. This correction was calculated first using the 
proximity force approximation 1651 . This approximation assumes that the surface 
profile slowly varied in comparison with the distance between the bodies. The lat- 
eral size of a rough profile is given by the correlation length ^ , therefore, PFA can 
be applied if £, ^ d. This condition is very restrictive since typical values of £, for 
deposited metals films (grain size) are in the range 20-100 nm. In most of the exper- 
imental situations the condition £, ^ d is broken. It was realized for the first time 
in Ref. iflTl . More general theory iflTl [Tsl [191 for the roughness correction can be 
applied at ^ < d and treats the correction perturbatively within the scattering for- 
malism (see Lambrecht paper in this volume). Here we discuss application of this 
theory to realistic rough surfaces and describe situations, for which one has to go 
beyond the perturbation theory to find agreement with experiments. 



3.2.1 Application of the perturbation tiieory 

Let us consider two parallel rough plates. A plate surface can be described by the 
roughness profile lii{x,y) {i =1,2 for plate 1 or 2) as shown in Fig. |9|a). The av- 
eraged value over large area is assumed to be zero {hi{x,y)) — 0. Then the local 
distance between the plates is 

d{x,y) ^ d - hi{x,y) - h2{x,y). (25) 

This distance depends on the combined rough profile h{x,y) = hi{x,y) + h2{x,y). 
As explained in Sec. 13.31 the interaction of two rough plates is equivalent to the 
interaction of a smooth plate and a rough plate with the roughness given by the 
combined profile li{x,y) (see Fig.|9|b)). 

Let us assume further that the interaction energy per unit area of two flat plates 
is Epp{d). If the rms roughness of the combined profile li{x,y) is small, w ^d, but 
the correlation length is large, ^ ^ J, we can present the interaction between rough 
plates as 

£;7" = (Epp (di^M « Epp{d) + ^ (h^) , (26) 

where = = Wj + Wj. Equation ( |27] i defines the PFA roughness correction 
5Epp = EppW^/2. This correction was used in all early studies to estimate the rough- 
ness effect. 

It was noted ifTTl that in most experimental configurations the condition ^ ^ d is 
broken and PFA cannot be applied. In Refs. lfT8l[T9l was developed a theory, which 
is not restricted by the condition ^ ^ t/. Within this theory the roughness correction 
is expressed via the spectral density of the rough surface G{k) as 

5Epp{d)= j -^G{k,d)a{k), (27) 
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Fig. 9 Contact of two rough surfaces, (a) Two rough plates in contact. Roughness of each plate, 
hi{x,y), is counted from the mean plane shown by the dashed line. The distance between bodies 
is the distance between these mean planes, (b) The interaction between two rough plates is equiv- 
alent to the interaction between a smooth plate and a rough plate with the roughness given by the 
combined profile h{x,y). The distance upon contact, da has well defined meaning in this case. See 
Sec. Upland Ref. ffS for details. 



where G{k,d) is a roughness response function derived in [T9^|. The PFA re- 
sult (IZTI i is recovered from here in the Hmit of small wavenumbers k ^ when 
G{k, d) — ^ Epp{d) /2. The roughness power spectrum is normalized by the condition 
/ d^k(y{k) / {In)^ ~ w^. The spectrum itself can be obtained from AFM scans and 
in the case of self-affine rough surfaces is given by Eq. (l24l i. 

Let us enumerate the conditions at which Eq. ( l27b is valid, (i) The lateral di- 
mensions of the roughness ^ must be much smaller than the system size L, ^ <C L. 
This is usually the case in experiments, (ii) The rms roughness w must be small 
compared to the separation distance, w <C (i. This condition means that roughness is 
treated as perturbative effect, (iii) The lateral roughness dimensions must be much 
larger than the vertical dimensions, w <C |19|. The last two assumptions are not 
always satisfied in the experiment. 

In the plate-plate configuration the force per unit area can be calculated as 
the derivative of E'pp^'\d). For the sphere-plate configuration, which is used in 
most of the experiments, the force is calculated with the help of PFA as Fsp{d) = 
2nREp'p^''{d). In contrast with the roughness correction the latter relation is justi- 
fied for d R, which holds true for most of the experimental configurations. We 
use the sphere-plate configuration to illustrate the roughness effect. The deposited 
gold films can be considered as self-affine. For all calculations reported here we 
are using our smoothest spheres with the parameters vv = 1 .8 nm, £, ^22 nm, and 
H = 0.9. We alter only the plate roughness since it is easy to prepare and replace 
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during experiments. We use the optical data for gold films described previously. It 
was found that the PFA limit is quickly recovered for increasing correlation length. 
Deviations from PFA prediction for real films were found to be about 1-5% in the 
range li = 50 — 200 nm. 

Therefore, for real rough surfaces the scattering theory gives a few percent cor- 
rection to the force compared to the PFA. This difference is difficult to measure. 
However, at small separations both PFA and perturbation theory fail since the rms 
roughness becomes comparable in size to the separation distance. It would be in- 
teresting to calculate the roughness effect when d is comparable with w. At the 
moment there is no a theoretical approach to estimate the effect except a direct nu- 
merical analysis similar to that used in 16611 . It will therefore be interesting to do a 
full numerical analysis for films with high local slopes instead of using perturbation 
theory. On the other hand it is experimentally possible to go to sufficiently small 
distances as it will be discussed below. 



3.2.2 Experimental evidence of large roughness effect 

The Casimir forces between a 100 /im gold coated sphere and substrates cov- 
ered with Au to different thicknesses from 100 nm to 1600 nm were measured in 
ll24ll . Different layers of Au resulted to different roughnesses and different correla- 
tion lengths, which are collected in Tab. [3] The roughness exponent was constant 
H = 0.9 ± 0.05 in agreement with former growth studies of thin films 16411621 . The 
sphere was attached to a cantilever with a spring constant of 0.2 N/m. The calibra- 
tion procedure is described in |24 1 . 

The force results are shown in Fig.[TO] Our measurements were restricted to sep- 
arations below 200 nm where the Casimir force is large enough compared to the 
approximately linear signal from laser light surface backscattering l26l l24ll . The 
small separation Umit or contact point is restricted by the jump to a contact 5 
nm) I67ll68ll69l and surface roughness. Note that an error of 1.0 nm in absolute dis- 
tance leads to errors in the forces of up to 20% at close separations as for example 
J ^ 10 nm fS?]. Thus we cannot detect the scattering effects described above. What 
we do see is the failure of the perturbation theory, for the roughest films, for which 
the roughness strongly increases the force. These deviations are quite large, result- 
ing in much stronger forces at the small separations < 70 nm. At larger separations, 
70-130 nm (within our measurement range), where the roughness influence is neg- 
ligible, the usual 1 /cf-^ scaling of the force observed also in other experiments with 
gold is recovered and agreement with the theory is restored. For the smoother films 
deviations from theory below 40 nm can be explained with the error in the distance. 

Qualitatively the roughness effect can be reproduced by calculating the force 
between small areas on the surfaces separated by the local distance d{x,y). One can 
call this procedure the non-perturbative PFA approach. Although it is qualitative, 
it can be used to obtain an estimate of the force at close proximity (2 nm above 
the point upon contact), where the roughness has an enormous influence on the 
Casimir force (see inset in Fig.lToli. This explains the jump to contact only partially. 
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Fig. 10 Casimir forces measured for the films of different roughness. The roughness effect man- 
ifests itself as a strong change in scaling at smaller separations, where the forces become much 
stronger. Errors in separation are shown for some points by the horizontal bars. The theoretical 
curves are shown for the 100 nm (smooth) and 1600 nm (rough) films. The inset shows the forces 
calculated by integrating over the roughness scans using PFA (see text). 



since approximately 5 nm above the point of contact, the capillary force will act as 
well. This force appears due to absorbed water and Kelvin condensation ITOl l22l 
resulting in water bridges formation between bodies. In the limit of fully wetted 
surface (see Fig.fTTli the capillary force is given by Fcap ~ AnjRco?, & (upper dashed 
line), while for a single asperity (of size £, ) wetting the minimum capillary force is 
Fcap ~ 4nY^ cos (lower dashed line). Here y is the surface tension of liquid, and 
I? is the contact angle 16911671 . 

At this point we have to compare the Casimir adhesion between rough films with 
adhesion by capillary forces Il68ll69l . While Casimir forces may lead to stiction be- 
tween movable parts, once the surfaces are in contact capillary forces (being present 
in air between hydrophilic surfaces into close proximity) are much stronger. The 
roughness effect on capillary adhesion is also much stronger as shown in Fig. [TT] 
Note that the Casimir force for a /? = 50 fim sphere is in the order of 10 nN at 10 nm 
separations. Capillary forces between a mica substrate and the same sphere are as 
large as 10 jxN deeming contact measurements with soft cantilevers (spring constant 
< 1 N/m) even impossible since the retraction range is outside of that of most piezo 
z-ranges. Measurements of the capillary forces with a smooth sphere used for the 
Casimir force measurement are shown in Fig.IH] One can see that when roughness 
increases a few times the force decreases by more than a factor of 100. This can be 
related to full sphere wetting and asperity wetting. The size of the sphere R — 50 jim 
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Fig. 11 Capillary forces in air (relative humidity 2-60%) for a smooth gold coated sphere, vv,p;, ~ 1 
nm, measured with a stiff cantilever, k = 4 N/m, and different rough films. The inset shows a fully 
wetted sphere (upper dashed line), and a roughness asperity wetted sphere (lower dashed line). 



is 1000 times larger than that of an asperity ^ 50 nm. Multiple asperity capillary 
bridge formation is likely to happen in the rough regime giving increasing forces. 

Furthermore, formation of capillary bridges means that under ambient conditions 
gold absorbs water, and as a result it is covered with an ultra thin water layer The 
experiments 1691 l67l suggest that the thickness of this layer is in the nanometer 
range, 1-2 nm. The natural questions one could ask is how thick the water layer is, 
and what is the influence of this water layer on the dispersive force ItTI ? At short 
separations, < 20 nm, these questions become of crucial importance because they 
place doubts in our understanding of the dispersive forces when experiments under 
ambient conditions are compared with predictions of the Lifshitz theory. Figure [12] 
shows the Casimir force measured at short distances together with theoretical calcu- 
lations made for gold with or without a water layer on top. The errors are shown to 
arise mainly from the experimental uncertainty in determining the separation upon 
contact do due to nanoscale surface roughness. We can conclude that the experiment 
can exclude the water layer thicker than 1.5 nm. Figure fTSl b) shows that the effect of 
water becomes very significant at separations below 10 nm, which were not acces- 
sible in our measurements due to jump-into-contact. We presented only the forces 
between flat surfaces because at these small separations there is no a reliable way to 
estimate the roughness correction. 
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Fig. 12 (a) Experimental data for the force vs distance (circles) down to 13 nm separations, and the 
theoretical prediction without water layer (red curve). Errors in the absolute separation are shown 
for some points by the bars. The continuous black curve is the prediction for continuous water 
layer of thickness h = 1.5 nm. The dashed black curve corresponds to the same water layer with 
50% of voids (/ = 0.5). (b) Theoretical calculations using the Lifsthiz theory for flat surfaces and 
continuous water layer films with thicknesses h =0, 1, 1.5 nm for small separations. 



3.3 Distance upon contact 

The absolute distance separating two bodies is a parameter of principal importance 
for the determination of dispersive forces. It becomes difficult to determine when 
the separation gap approaches nanometer dimensions. This complication originates 
from the presence of surface roughness, which manifests itself on the same scale. In 
fact, when the bodies are brought into gentle contact they are still separated by some 
distance do, which we call the distance upon contact due to surface roughness. This 
distance has a special significance for weak adhesion, which is mainly due to van der 
Waals forces across an extensive noncontact area 1211 . It is important for MEMS and 
NEMS as unremovable reason for stiction ll22l . In the modern precise measurements 
of the dispersive forces do is the main source of errors (see reviews ||5]|6l). This 
parameter is typically determined using electrostatic calibration. The distance upon 
contact is usually considerably larger than the rms roughness because it is defined 
by the highest asperities. It is important to clear understand the origin of do, its 
dependence on the lateral size L of involved surfaces, and possible uncertainties in 
its value 1721 . These are the questions we address in this subsection. 



3.3.1 Plate-plate contact 

Two plates separated by the distance d and having roughness profiles hi{x,y) locally 
are separated by the distance d{x,y) given by Eq. (|25] | (see Fig. |9]l. Indeed, the 
averaged local distance has to give d, {d{x,y)) = d. We can define the distance upon 
contact do as the largest distance d ~ do, for which d{x,y) becomes zero. 



26 



P. J. van Zwol, V.B. Svetovoy, and G. Palasantzas 



It is well known from contact mechanics 1731 that the contact of two elastic rough 
plates is equivalent to the contact of a rough hard plate and an elastic flat plate with 
an effective Young's modulus E and a Poisson ratio v. Here we analyze the contact 
in the limit of zero load when both bodies can be considered as hard. This limit 
is realized when only weak adhesion is possible, for which the dispersive forces 
are responsible. Strong adhesion due to chemical bonding or due to capillary forces 
is not considered here. This is not a principal restriction, but the case of strong 
adhesion has to be analyzed separately. Equation dZSl l shows that the profile of the 
effective rough body is given by 



The latter means that h{x,y) is given by the combined image of the surfaces facing 
each other. If topography of the surfaces was determined with AFM, we have to 
take the sum of these two images and the combined image will have the size of the 
smallest image. 

To determine c/q we collected 1721 high resolution megascans (size 40 x 40 /im^, 
lateral resolution 4096 x 4096 pixels) for gold films of different thicknesses de- 
scribed before. The maximal area, which we have been able to scan on the sphere, 
was 8 X 8 lirn^ (2048 x 2048 pixels). The images of 100 nm film, sphere, and 1600 
nm film are shown in Fig.[T|(a), (b), and (c), respectively. Combining two images 
and calculating from them the maximal peak height we can find do for a given size 
of the combined image. Of course, taking the images of the same size every time we 
will get different value of dd and averaging over a large number of images we will 
find the averaged do and possible rms deviations. This is quite obvious. What is less 
obvious is that if we take images of different size and will do the same procedure 
the result for c/q will be different. 

Let Lq be the size of the combined image. Then, in order to obtain information 
on the scale L = Lq/2", we divide this image on 2" subimages. For each subimage 
we find the highest point of the profile (local c/q), and average all these values. This 
procedure gives us doiL) and the corresponding statistical error. Megascans are very 
convenient for this purpose otherwise one has to collect many scans in different 
locations. For the 100 nm film above the 400 nm film the result of this procedure 
is shown in Fig. [T3l a). We took the maximum area to be 10 x 10 /im^. The figure 
clearly demonstrates the dependence of do on the scale L although the errors appear 
to be significant. The dependence of the rms roughness w on the length scale L is 
absent in accordance with the expectations, while only the error bars increase when 
L is decreasing. 

To understand the dependence do (L) let us assume that the size L of the area of 
nominal contact is large in comparison with the correlation length, L ^ ^ . It means 
that this area can be divided into a large number A^^ = L?' /^^ of cells. The height 
of each cell (asperity) can be considered as a random variable h. The probability to 
find h smaller than some value z can be presented in a general form 



h{x,y)=hi{x,y)+h2{x,y). 



(28) 



(29) 
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Fig. 13 (a) Distance upon contact as a function of the length scale. Dots with the error bars are 
the values calculated from the megascans. The solid curve is the theoretical expectation according 
to Eq. ([30). Note that is considerably larger than w. (b) Statistics of the surface roughness. 
Four 10 X 10 llm^ images were used. The main graph shows the "phase" as a function of z- The 
green lines show the best fits at large positive and negative z- The dashed red lines demonstrate the 
solution of Eq. <30t . The inset shows the cumulative distribution P{z). 



where the "phase" (z) is a nonnegative and nondecreasing function of z- Note that 
( [29] ) is just a convenient way to represent the data: instead of cumulative distribu- 
tions P{z) we are using the phase (/>(z), which is not so sharp function of z- 

For a given asperity the probability to find its height above c/q is 1 —Pido), then 
within the area of nominal contact one asperity will be higher than do if 



lV^') = 1 or .^)(^o)=ln(LV^')- (30) 



This condition can be considered as an equation for the asperity height because due 
to a sharp exponential behavior the height is approximately equal to do. To solve 
( [30l l we have to know the function 0(z), which can be found from the roughness 
profile. 

The cumulative distribution P{z) can be extracted from combined images by 
counting pixels with the height below z- Then the "phase" is calculated as (j){z) = 
— ln(l — P). The results are presented in Fig. [TSlb). The procedure of solving Eq. 
([30b is shown schematically in Fig.[T3jb) by dashed red curves, and the solution it- 
self is the red curve in Fig. [T3l' a). It has to be mentioned that the normal distribution 
fails to describe the data at large z- Other known distributions cannot satisfactory 
describe the data for all z- Asymptotically at large |z| the data can be reasonably 
well fit with the generahzed extreme value Gumbel distributions (green lines in Fig. 

[Hb)) m-. 

The observed dependence do{L) can be understood intuitively. The probability to 
have one high asperity is exponentially small but the number of asperities increases 
with the area of nominal contact. Therefore, the larger the contact area, the higher 
probability to find a high feature within this area. Scale dependence of do shows 
that smaller areas getting into contact will be bound more strongly than larger areas 
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because upon the contact they will be separated by smaller distances. This fact is 
important for weak adhesion analysis. 

3.3.2 Sphere-plate contact 

Most of the Casimir force experiments measure the force in the sphere-plate con- 
figuration to avoid the problem with the plates parallelism. Let us consider how the 
scale dependence of do manifests itself in this case. Assuming that the sphere is 
large, d, the local distance is 

d{x,y) ^d+{x'- +/) /2R - h{x,y), (32) 

where h{x,y) is the combined topography of the sphere and the plate. As in the 
plate-plate case do is the maximal d, for which the local distance becomes zero. 
This definition gives 

t/o=max[/z(x,y)-(x2+/)/27?l . (33) 

.V,V V / J 

Now do is a function of the sphere radius R, but, of course, one can define the length 
scale Lr corresponding to this radius R. 

As input data in Eq. (l33T l we used the combined images of the sphere and different 
plates. The origin {x = 0, y = 0) was chosen randomly in different positions and then 
do was calculated according to (l33T l. We averaged do found in 80 different locations 
to get the values of d'Q', which are collected in Tab. |3] The same values can be 
determined theoretically using do{L) found between two plates (see Eq. (l30l l and 
Ref. (721). 





100 nm 


200 nm 


400 nm 


800 nm 


1600 nm 


w 
d'o' 


3.8 

26.1(3.8) 
12.8(2.2) 
17.7(1.1) 


4.2 

28.8(3.7) 
15.9(2.7) 
20.2(1.2) 


6.0 

34.4(4.7) 
24.5(4.8) 
23.0(0.9) 


7.5 

30.6(2.4) 
31.3(5.4) 
34.5(1.7) 


10.1 

42.0(5.5) 
55.7(9.3) 
50.8(1.3) 



Table 3 The pai'ameters characterizing the sphere-film systems (all in nm) for S = 50 /im. The 
first three rows were determined from combined images. The last row for gives the values of 
do determined electrostatically. The errors are indicared in brackets. 



The values of c/™ for rougher films are in agreement with those found from the 
electrostatic calibration. However, for smoother films (100 and 200 nm) d'^' and dQ 
do not agree with each other. This is most likely to be attributed to the roughness 
on the sphere, which varies considerably locally. For example, between those 80 
found in different locations 5% are in agreement with c/q' found from the elec- 
trostatic calibration 1721 . This is illustrated by the fact that when the roughness of 
the plate dominates the discrepancy between t/™ and d^ disappears. Note that the 
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Fig. 14 Schematic explanation of additional uncertainty 8d in do (see text). The sphere in two 
positions is shown by the dashed (contact) and solid blue curves. The inset shows the variance of 
dd as a function of the scale L or separation d. 



Standard deviations for t/™ are larger than that for . The standard deviations in 
d™ originate from place to place variations of c/q". In the case of electrostatic deter- 
mination of diy statistical variation of c/q from place to place is not included in the 
errors of d^ . This explains why the errors in d^ are smaller 

Consider the experimental situation when the dispersive force is measured in 
the sphere-plate configuration. The system under consideration is equivalent to a 
smooth sphere above a combined rough profile h{x,y). The position of the average 
plane depends on the area of averaging especially for small scales L. The profile 
shown in Fig. [14] demonstrates different mean values in the left and right segments 
shown by the dashed black lines. Both of these values deviate from the middle line 
for the scale 2L (solid black line). The true average plane is defined for L — oo. 

Position of the average plane define the absolute separation of the bodies. It has 
to be stressed that the electrostatic and Casimir interactions "see" different areas 
on the plate. This is due to different dependence on d and quite often the electro- 
static calibration is performed at larger separations than measurement of the Casimir 
force. The size L of the interaction area is determined by the relation = anRd, 
where a = 2 for the electrostatic and a = 2/3 for pure Casimir force. Therefore, 
these two interactions can "see" different positions of the average planes. It intro- 
duces an additional uncertainty 8d in the absolute separation ll72l . For a fixed L 
this uncertainty is a random variable distributed roughly normally around 5d = 0. 
However, it has to be stressed that dd manifests itself not as a statistical error but 
rather as a kind of a systematic error This is because at a given lateral position of 
the sphere this uncertainty takes a fixed value. The variance of 5d is defined by the 
roughness statistics. It was calculated from the images and shown as inset in Fig. [14] 
One has to remember that with a probability of 30% the value of 5d can be larger 
than that shown in Fig. [14] 
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In this chapter we considered the Casimir force between realistic materials con- 
taining defects, which influence the optical properties of interacting materials, and 
having surface roughness, which contributes to the force. 

It was demonstrated that the gold films prepared in different conditions have 
dielectric functions, which differ from sample to sample, and this difference cannot 
be ignored in the calculation of the Casimir force aimed at precision better than 10%. 
The main conclusion is that for metals one has to measure the dielectric function of 
used materials in a wide range of frequencies, where far and mid IR are especially 
important. Precise knowledge of the dielectric functions is also important for low 
dielectric materials. In this case significant sensitivity of the force to the dielectric 
functions is realized nearby the attractive-to-repulsive transition in solid-liquid-solid 
systems. 

The roughness correction to the Casimir force can be reliably calculated if rms 
roughness w is small in comparison with the separation, w <C t/, when one can apply 
the perturbation theory. In the experiments at short separations this condition can be 
violated. The current situation with the theory is that there is no direct method to cal- 
culate the force between rough bodies when d except using rather complicated 
numerical calculations. 

We gave also a detailed analysis of the distance upon contact, cIq, which is an 
important parameter in Casimir physics. Analysis of AFM scans demonstrated that 
do is always a few times larger than the rms roughness. Moreover, c/q is a function 
of the size L of the nominal area of contact. This dependence is important for weak 
adhesion, which is due to van der Waals forces across an extensive noncontact area. 
Uncertainty in do is the main source of errors in the Casimir force measurements. 
We demonstrated here that there is an additional indefiniteness in do, which cannot 
be excluded by the electrostatic calibration. It becomes very important for small 
areas of interaction. Also, this indefiniteness has to be taken into account if one 
compares two independent experiments. 
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